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Abstract. The paper presents an investigation of aerosol 
seasonal variations in several urban-industrial regions in 
Ukraine. Our analysis of seasonal variations of optical and 
physical aerosol parameters is based on the sun-photometer 
2008-2013 data from two urban ground-based AERONET 
(AErosol RObotic NETwork) sites in Ukraine (Kyiv, Lu- 
gansk) as well as on satellite POLDER instrument data for 
urban-industrial areas in Ukraine. We also analyzed the data 
from one AERONET site in Belarus (Minsk) in order to 
compare with the Ukrainian sites. Aerosol amount and opti- 
cal depth (AOD) values in the atmosphere columns over the 
large urbanized areas like Kyiv and Minsk have maximum 
values in the spring (April-May) and late summer (August), 
whereas minimum values are observed in late autumn. The 
results show that fine-mode particles are most frequently de- 
tected during the spring and late summer seasons. The analy- 
sis of the seasonal AOD variations over the urban-industrial 
areas in the eastern and central parts of Ukraine according 
to both ground-based and POLDER data exhibits the similar 
traits. The seasonal variation similarity in the regions denotes 
the resemblance in basic aerosol sources that are closely re- 
lated to properties of aerosol particles. The behavior of basic 
aerosol parameters in the western part of Ukraine is differ- 
ent from eastern and central regions and shows an earlier ap- 
pearance of the spring and summer AOD maxima. Spectral 
single-scattering albedo, complex refractive index and size 
distribution of aerosol particles in the atmosphere column 
over Kyiv have different behavior for warm (April-October) 


and cold seasons. The seasonal features of fine and coarse 
aerosol particle behavior over the Kyiv site were analyzed. A 
prevailing influence of the fine-mode particles on the optical 
properties of the aerosol layer over the region has been estab- 
lished. The back-trajectory and cluster analysis techniques 
were applied to study the seasonal back trajectories and pre- 
vailing directions of the arrived air mass for the Kyiv and 
Minsk sites. 


1 Introduction 

Aerosol seasonal variations have been investigated in differ- 

ent regions with various ground-based (Gerasopoulos et al., 
2007; Jaroslawski and Pietruczuk, 2010; Rana et al., 2009; 
Andrews et al., 2011; Leskinen et al., 2012; Liu et al., 2012; 
Pietruczuk and Chaikovsky, 2012), satellite (e.g., Barnaba 

and Gobbi, 2004; Song et al., 2008), combined ground- 
based and satellite (e.g., Zawadzka et al., 2013; Bovchaliuk 
et al., 2013), and balloon-borne (Hara et al., 2011, 2013) 
techniques for aerosol measurements. The NASA AErosol 
RObotic NETwork (AERONET; Holben et al., 1998; http: 
//aeronet.gsfc.nasa.gov/) data are used intensively in such 

studies. For example, AERONET sun-photometer observa- 
tions have recently been used (Liu et al., 2012) for study 
of the seasonal variations in aerosol optical properties in 
China, including the aerosol optical depth (AOD), Angstrom 
exponent (AE), and single-scattering albedo (SSA). For the 
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eastern China region it has been concluded that the AOD is 
largest in the summer and smallest in the winter, whereas 
the SSA values exhibit weak seasonal variations with the 
smallest values during the winter and the largest during the 
summer. The seasonal behavior of aerosol optical properties 
determined from vertically resolved in situ measurements 
over rural Oklahoma, USA, were compared with AERONET 
data by Andrews et al. (2011). Combined ground-based 
(AERONET/PHOTONS) and satellite (Moderate-Resolution 
Imaging Spectrometer, or MODIS) data have been used 
to study seasonal aerosol content and properties over Eu- 
rope for 2000-2009 and their impact on ultraviolet (UV) 
radiance (Chubarova, 2009). The analysis by Jaroslawski 
and Pietruczuk (2010) of seasonal AOD variations at Belsk 
(Poland) based on data over the 2002-2007 period reveals 
two AOD maxima: in April and July-August. In the paper 
by Zawadzka et al. (2013), the difference in AOD between an 
urban area (Warsaw) and a rural area (Belsk) during 2005- 
2011 was investigated. The estimated influence of urban-area 
emissions on the AOD was reported as being less than 15 %. 
As reported by Uscka-Kowalkowska (2013), the extinction 
of direct solar radiation has been studied during the period 
1964-2003 in an elevated mountainous region of Poland that 
can be considered as a non-urban setting. The increased ex- 
tinction was observed during the spring and summer seasons 
and showed traits similar to those in urban-industrial areas of 
eastern Europe. This suggests that the behavior of the aerosol 
content in this region is governed mainly by the aerosol trans- 
port. 

In previous analyses (e.g., see Chubarova, 2009; Giles et 
al., 2012; Pietruczuk and Chaikovsky, 2012; Zawadzka et 
al., 2013), the eastern European region is considered to be 
the source of urban-industrial aerosols according to the gen- 
eral aerosol type classification by Dubovik et al. (2002). In- 
deed, there are many existing and potential aerosol pollution 
sources in this region, including Ukraine: intensive trans- 
port and agriculture, heavy industry, open steppe fields, and 
mining. Furthermore, this region is characterized by numer- 
ous forest, grassland, and peat wildfires, and also can some- 
times be considered as a source of biomass burning aerosols 
(Bamaba et al., 2011; Witte et al., 2011; Bovchaliuk et al., 
2013). The steppe regions of southern Ukraine also experi- 
ence periodical dust storms (e.g., Birmili et al., 2008). 

The investigated region is mostly fiat terrain with altitude 
difference between the center and east of Ukraine of about 
300 m above sea level (a.s.l.). The developed heavy indus- 
try is concentrated in the Donetsk, Lugansk, and Kharkiv 
metropolitan areas; industry and open mines are typical of 
east-central Ukraine (Dnipropetrovsk and Kirovograd areas); 
while steppes are located in the south and the southeast of 
Ukraine (Fig. 1). 

In the west of Ukraine the altitude increases to about 
400 m a.s.l. with mountains up to 2000 m high in the far 
southwest. This part of Ukraine is less populated and 
has lower industrial load (the Rivne and Lviv city areas). 
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Figure 1. Locations of the AERONET sites (Kyiv, Lugansk, Sev- 
astopol, and Minsk) and areas covered by POLDER measurements 
(Kyiv, Lugansk, Donetsk, Kharkiv, Dnipropetrovsk, Lviv, Rivne, 
and Minsk) as discussed in this paper. 

Therefore, the potential urban-industrial aerosol sources are 
concentrated in the east, east-central and southeast regions. 
For big cities like Kyiv and Kharkov, a large potential source 
of aerosols is car traffic. For example, the concentration of 
cars in the central part of Kyiv sometimes reaches 800 vehi- 
cles per kilometer, producing ~ 90 % of all atmospheric pol- 
lution in the area. 

Previous studies (e.g., Chubarova, 2009; Bovchaliuk et 
al., 2013; Ogunjobi et al., 2008; Israelevich et al., 2012) 
show that the aerosol seasonal behavior in the eastern Eu- 
ropean region strongly depends on regional/latitudinal con- 
ditions. The seasonal variability of aerosols in the atmo- 
sphere over Ukraine has been less studied due to the lack 
of data collection from ground-based observations in previ- 
ous years. A preliminary analysis of general aerosol param- 
eters for the eastern European region has been done by the 
authors (Bovchaliuk et al., 2013). In continuation of our pre- 
vious study, we focus here on an analysis of seasonal aerosol 
properties and dynamics over the urban-industrial areas of 
Ukraine based on AERONET sites and POLDER/PARASOL 
satellite imaging spectroradiometer data. Unlike the preced- 
ing paper (Bovchaliuk et al., 2013) we analyze the variation 
of the total AOD, the fine- and coarse-mode components of 
the AOD and volume concentration, the spectral SSA, and 
the refractive index (RI) for AERONET sites in Ukraine and 
expand the analysis of the POLDER data to a larger number 
of urban areas. The purpose of the analysis is to determine 
dominant aerosol types in individual seasons. In Sect. 2, we 
describe the instrumentation and measurement sites/regions. 
Seasonal variations of aerosol parameters are analyzed in 
Sect. 3. We discuss and compare the results for individual 
sites and regions in Sect. 4 and make concluding remarks in 
Sect. 5. 
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2 Description of instrumentation and 
measurement sites 

2.1 AERONET data and sites 

A convenient and powerful method for studying the aerosol 
content and parameters is ground-based measurements 
with the AERONET/PHOTONS sun-photometer network 
(Holben et al., 1998), which allows obtaining long-term se- 
ries of uniform and accurate data for the analysis of varia- 
tions on different timescales. The sun-photometer data are 
also useful in satellite retrieval validation (e.g., see King et 
al., 1999; Fan et al., 2008; Li et al., 2009; Mishchenko et 
al., 2010; Schuster et al., 2012). Studies of aerosol content 
and property data collected by satellite measurements over 
specific regions with AERONET sites and joint satellite and 
ground-based data analyses allow one to determine aerosol 
seasonal dynamics and deduce the aerosol behavior pattern. 
We use this approach in our investigations for the Ukrainian 
territory. The AERONET network has been recently ex- 
panded in the central part of eastern Europe (Danylevsky 
et al., 2011; Milinevsky et al., 2012; Milinevsky, 2013). At 
present, the two AERONET sites are continuously operated 
there: the Kyiv site (operated since 2008) in Ukraine and the 
Minsk site (operated since 2002) in Belarus. The data from 
the Sevastopol site (operated since 2006) in Ukraine have 
been described in detail by Bovchaliuk et al. (2013) and are 
partly used in our analysis in Sect. 4. Also data from the Lu- 
gansk site, which has been operated periodically since 2011, 
are included in our consideration in comparison to POLDER 
data. The data from these three AERONET sites in Ukraine 
and one site in Belarus for comparison purposes have been 
used for analysis of aerosol seasonal features in this paper 
together with POLDER data for several urban-industrial re- 
gions in Ukraine. The information on AERONET sites and 
the Ukrainian urban-industrial areas selected for our anal- 
ysis using POLDER data is given in Table 1. The Kyiv 
AERONET site has the longest data record in the region. 
Therefore, to evaluate seasonal aerosol variations we perform 
the most detailed analysis for that area. 

For the purposes of this paper, mostly level 2.0 cloud- 
screened and quality-assured data were used. In a few cases, 
when level 2.0 data were not available, we used level 1.5 
(cloud-screened and outliers removed) data. In such cases the 
sign “level 1.5” appears in relevant figure. The AERONET 
Version 2 direct sun and inversion algorithms and software 
were used for the determination of aerosol properties from 
sun-photometer measurements (Dubovik and King, 2000; 
Dubovik et al., 2000, 2006; http://aeronet.gsfc.nasa.gov/). 

For aerosol dynamics analysis in this paper, we extracted 
the aerosol Version 2 direct sun algorithm (DSA) level 
2.0 parameters from the AERONET database: daily and 
monthly weighted averages of AOD at 440 and 870 nm, 
and the Angstrom exponent values computed using the 440- 
675 nm and 440-870 nm channels. The almucantar retrievals 


Inversion Products level 2.0 and level 1.5 were obtained from 
the AERONET database: level 2.0 size distribution and vol- 
ume concentration, level 1.5 SSA and complex RI, level 1.5 
absorption AOD and extinction AOD. To determine the sea- 
sonal behavior, we calculated monthly averaged parameters 
during the period 2008-2013 for the Kyiv site, 2002-2013 
for the Minsk site, and November 2011-August 2013 for 
the Lugansk site. The Angstrom exponent (675-870 nm) was 
calculated from level 2.0 AOD data. 

The AOD value is a general measure of the total aerosol 
amount in the atmospheric column. The AE is a parame- 
ter related to the aerosol particle size distribution and de- 
pends on the effective radius of particles that allows divid- 
ing aerosol by fraction: a small AE pertains to coarse-mode 
particles, while a large AE corresponds to small (fine-mode) 
particles. But it is influenced also by effective variance of the 
size distribution (Tanre et al., 2001; Schuster et al., 2006; 
Kokhanovsky, 2008). For the AE value analysis, we cal- 
culated the AE according to the AERONET procedures as 
a(ki, k 2 ) =ln[AOD(Ai)/AOD(A 2 )]/ln(A 2 Ai) using AOD 
data (Holben et al., 1998). 

We use also the volume concentration AERONET pa- 
rameter, which is the total volume of aerosol particles per 
unit surface area (in pm 3 pm -2 ), to infer its seasonal vari- 
ation. The volume concentration value is calculated from 
sun-photometer measurements of the spectral sky radiance 
distribution along the solar almucantar using the inversion 
algorithm and software by Dubovik and King (2000). The 
standard AERONET inversion algorithm provides these data 
under the assumption of a bimodal log-normal particle size 
distribution using sky radiance measurements along the solar 
almucantar for solar zenith angles (SZA) ranging from 75 to 
50°. 

The accuracy of retrieving aerosol particle properties is 
an important aspect of solving inverse problems. The spec- 
tral AOD values from the AERONET Version 2 DSA are 
retrieved with errors ±0.01 in the visible and near-infrared 
parts of the spectrum and with a larger uncertainty (±0.02) in 
the UV band (Holben et al., 1998). The AE uncertainty is on 
the order of ~ 0. 1 . The retrieval accuracy for the particle vol- 
ume size distribution dV(r)/dlnr, complex RI, and SSA de- 
pends on many factors, e.g., the sun-photometer calibration, 
AOD, solar zenith angle range, and particle sizes (Dubovik et 
al., 2000). The accuracy of the particle volume size distribu- 
tion dV (r) /din r is adequate essentially in the entire range of 
AODs (e.g., AOD (440 nm) > 0.05). In particular, for the in- 
termediate particle size range (0.1 < r < 7 pm), the retrieval 
errors do not exceed 10% in the maxima of the size distri- 
bution and may increase up to about 35 % for the size bins 
corresponding to the minimum values of dV (r) /din r in this 
size range. The dV (r)/dlnr retrieval errors rise to 80-100 % 
(and even higher) for sizes smaller than 0.1 pm and larger 
than 7 pm (Dubovik et al., 2000). 

The SSA values with an accuracy of ~ 0.03 and complex 
RI values with errors of ~ 30-50% for the imaginary part and 
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Table 1 . AERONET sites and Ukrainian urban-industrial areas selected for analysis. 


Site/area Coordinates, elevation a.s.l. Instrumentation 


Kyiv, Ukraine 
Lugansk, Ukraine 
Donetsk, Ukraine 
Kharkiv, Ukraine 
Dnipropetrovsk, Ukraine 
Rivne, Ukraine 
Lviv, Ukraine 
Sevastopol, Ukraine 
Minsk, Belarus 


50.36° N, 30.50° E, 200 m 
48.57° N, 39.36° E, 90 m 
48.03° N, 37.81° E, 300 m 
50.00° N, 36.23° E, 135 m 
48.45° N, 34.98° E, 68 m 
50.62° N, 26.25° E, 230 m 
49.83° N, 24.00° E, 270 m 
44.62° N, 33.52° E, 80m 
53.92° N, 27.60° E, 200 m 


AERONET/POLDER 

AERONET/POLDER 

POLDER-AERONET 

POLDER 

POLDER 

POLDER 

POLDER 

AERONET/POLDER 

AERONET/POLDER 


~ 0.04 for the real part can be retrieved only for high aerosol 
loading, with AOD (440 nm) >0.5. For observations with 
lower aerosol loadings, the retrieval accuracy of SSA and RI 
estimations significantly worsens. Dubovik et al. (2000) have 
shown that the accuracy levels worsen to 0.05-0.07 for SSA, 
80-100% for the imaginary part of the RI, and 0.05 for the 
real part of the RI when AOD (440 nm) < 0.2. Though the 
uncertainties in the aerosol characteristics over observational 
sites can be significant, we assume that monthly averaged 
aerosol parameters over the whole period of observation al- 
low us to describe the seasonal behavior of aerosol properties 
in the represented regions of Ukraine. 

2.2 POLDER/PARASOL data and regions 

Satellite instruments can provide the aerosol distribution 
data in the atmosphere over eastern Europe. The data ac- 
quired by the POLDER-3 imager on the PARASOL satellite 
(http://www.icare.univ-lillel.fr/parasol/; operated from De- 
cember 2004 to October 2013) have been used to investi- 
gate aerosol seasonal variations over the listed specific re- 
gions of Ukraine (Fig. 1 and Table 1). The aerosol parame- 
ters derived from the POLDER measurements are the AOD 
in the 865 nm band and the AE determined using the AODs in 
the 865 and 670 nm bands. The instrument allows obtaining 
the aerosol amount and the parameters of fine-mode parti- 
cles with sizes less than 0.3 pm from measurements of polar- 
ized back-scattered solar radiation (Deuze et al., 2001; Fan 
et al., 2008; Tanre et al., 2011). The fine-mode aerosols are 
produced mainly by anthropogenic sources and by biomass 
burning events, anthropogenic or natural (Tanre et al., 2001; 
Dubovik et al., 2002; Anderson et al., 2005; Kaufman et 
al., 2005; Amiridis et al., 2010; Chubarova et al., 2012). 
The polarized light measurements serve for improvement 
of the retrieval accuracy over land surfaces (Deuze et al., 
2001; Kokhanovsky et al., 2007; Mishchenko et al., 2007; 
Cairns et al., 2009), and this was the main reason for using 
POLDER instrument data for the analysis of the atmospheric 
contamination over Ukraine by anthropogenic aerosols. The 
POLDER instrument and data-processing algorithms are de- 
scribed by Deschamps et al. (1994), Breon (2006), and 


Kokhanovsky et al. (2007) (http://www.icare.univ-lillel.fr/ 
parasol/). The standard algorithm is based on the approach 
wherein the polarized light scattered by the aerosols is an- 
alyzed by assuming monomodal models of spherical non- 
absorbing particles and uses the look-up table method by 
Deuze et al. (2001). The retrieval accuracy has been in- 
vestigated for aerosol parameters over ocean (Herman et 
al., 2005; Goloub et al., 1999) and land. The correlation 
with AERONET sun-photometer data is 0.77-0.95 depend- 
ing on the location of ground-based sites (Breon et al., 2011; 
Fan et al., 2008; Bovchaliuk et al., 2013). The standard 
cloud-screening algorithm has been used for the processing 
of POLDER measurements over land surfaces (Breon and 
Colzy, 1999). In our investigation we use only cloud- free pix- 
els. More details on analysis methodologies using POLDER 
data can be found in our previous paper (Bovchaliuk et al., 
2013). 

The monthly mean AOD values at 865 nm and AE val- 
ues at 670-865 nm for the period from 2008 to 2011 over 
the Ukrainian urban-industrial areas around the cities of 
Kyiv, Lugansk, Donetsk, Kharkiv, Dnipropetrovsk, Lviv, and 
Rivne and the Belarus city of Minsk were derived from 
POLDER data using the algorithm described in Deuze et 
al. (2001). For the analysis, we used POLDER data for 
18 km x 18 km grid cells (Deuze et al., 2001; Bovchaliuk et 
al., 2013) on the Earth surface centered on cities mentioned 
above. 

It should be emphasized that satellite measurements are 
affected by the snow cover and significant cloudiness in 
the winter. Moreover, fewer data are available in the win- 
ter for latitudes above 50° N, because the solar elevation is 
too low to allow the AOD retrieval. This is the reason why 
the sparse data for January were not analyzed. The under- 
lying surface in this region varies in the spring and autumn 
due to the variable vegetation cover. In addition, the surface 
effect on polarization is relatively greater when the amount 
of aerosols in the atmosphere is small. However, for pur- 
poses of the POLDER and AERONET data comparison, it 
was considered that the environmental conditions at the Kyiv 
and Minsk AERONET sites are very similar. The sites are 
located nearby or within big cities surrounded by forests 


Atmos. Meas. Tech., 7, 1459-1474, 2014 


www.atmos-meas-tech.net/7/1459/2014/ 


G. Milinevsky et al.: Aerosol seasonal variations over urban-industrial regions 


1463 


and agricultural lands with similar, relatively flat landscapes. 
This allows us to assume that the reflectance properties of the 
land surfaces around these sites for POLDER retrievals are 
similar during a year. The comparison of the POLDER AOD 
(865 nm) and the AERONET AOD (870 nm) for fine-mode 
particles (r < 0.3 pm) measured at the Moscow, Minsk, Belsk, 
Kyiv, Moldova, and Sevastopol AERONET sites shows good 
agreement of the respective data. In general, the correla- 
tion coefficients for AOD (865nm)-AOD (870 nm) are in 
the range 0.76-0.93 except the value for the Sevastopol 
site (0.63), while the standard deviations in the POLDER- 
AERONET comparisons are less than 0.02 for the range of 
AOD values between 0.01 and 0.25. A detailed comparison 
of AERONET sun-photometer data and POLDER instrument 
data is given for all eastern European AERONET sites in 
Bovchaliuk et al. (2013). 


3 Results 

3.1 Seasonal variability of AOD and volume 
concentration 

For analysis of aerosol dynamics and optical properties in the 
atmosphere over eastern European region, the AOD (440 nm) 
AERONET sun-photometer data are most appropriate be- 
cause of the highest statistical accuracy of the measurements 
as compared with other spectral channels of the sun pho- 
tometer. In general, the data of the 440 nm spectral channel 
were often used for aerosol content and optical property stud- 
ies in different regions of the planet (e.g., Dubovik et al., 
2002; Boselli et al., 2012). Furthermore, the AOD (440 nm) 
is used as reference data in the AERONET algorithm in the 
accuracy estimations of the aerosol particle optical and mi- 
crophysical properties (Dubovik and King, 2000; Dubovik et 
al., 2000, 2006). 

However, it should be mentioned that the 440 nm spec- 
tral channel corresponds to the NO 2 molecular absorption 
band (Burrows et al., 1998). The NO 2 impact is consid- 
ered in the AERONET algorithms using the satellite data 
from OMI, SCIAMACHY, and GOME instruments (see 
AERONET Version 2 DSA, http://aeronet.gsfc.nasa.gov/). 
But the satellite NO 2 data uncertainties are still rather large 
(Boersma et al., 2004, 2011). Since the NO 2 concentration 
in the troposphere is determined mainly by anthropogenic 
sources like fossil fuel combustion (Enhalt et al., 2001), in a 
case of urban-industrial regions we therefore cannot exclude 
probable influence of NO 2 on measured AOD (440 nm) val- 
ues. To estimate this influence we can use the optical remote 
sensing data received for urban regions in USA (Alexandrov 
et al., 2002a, b). Though the yearly mean values of the NO 2 
atmospheric extinction, e.g., over New York city (measure- 
ments at 440 nm from September 1995 till November 1996) 
were below 5 Dobson units (DU), on some days the extinc- 
tion value reached 15 DU. Considering the absorption cross 


sections of NO 2 given in Burrows et al. (1998), the corre- 
sponding extinction optical thickness is 0.07 for the 5 DU 
and 0.2 for the 15 DU case. This value exceeds significantly 
the errors of AOD measurements by AERONET sun pho- 
tometers in the 440 nm spectral band (Holben et al., 1998). 

According Rublev et al. (2004) the NO 2 concentration at 
0.5 km altitude in Moscow at rush hour times corresponds to 
extinction of about 13.7 DU and that is approximately equal 
to the AOD (440 nm) value of about 0.18. The measurements 
over the Kyiv and Minsk sites show that mean values of the 
AOD (440 nm) are 0.25 and 0.23, respectively (see Table SI 
in the Supplement), and in most of cases (with few excep- 
tions of forest wildfire event and dust storm; Bovchaliuk et 
al., 2013) they are less than 0.5, which is a crucial value 
for the AERONET inversion algorithm to the SSA and RI 
aerosol particle determination. Though Kyiv and Minsk are 
not such big cities compared to New York and Moscow, 
NO 2 could impact AERONET AOD (440 nm) data in urban- 
industrial regions like those considered in our paper. Among 
the spectral channels of AERONET sun-photometer CE318- 
2 the channel 870 nm is free from NO 2 influence. Moreover 
this spectral band corresponds to the main spectral channel 
of the POLDER instrument. However, in Kyiv and Minsk the 
AOD (870 nm) are 2-3 times lower than the AOD (440 nm), 
and consequently the relative errors of the AOD (870 nm) 
measurements are proportionally larger than for the AOD 
(440 nm). Therefore, it looks reasonable to use the data from 
both spectral channels for evaluation of seasonal variations 
of aerosol optical properties in this paper. 

The monthly mean AOD (440 and 870 nm) values and AE 
(440-870 nm) data retrieved by the AERONET DSA have 
been averaged over the 2008-2013 period (the Kyiv site) and 
over the 2002-2013 period (Minsk site) to reveal the aerosol 
content and optical properties seasonal peculiarities. Accord- 
ing to our results, the AOD values exhibit two maxima: in 
spring (April-May) and in the second half of summer (Au- 
gust). Note that the summer AOD maximum is larger and is 
seen in both spectral channels, but spring maximum is clearly 
seen in 440 nm only and occurs mostly in April (the Kyiv 
site) and in April-May (the Minsk site) (Fig. 2a-c). 

The aerosol spectral characteristics in 440-870 nm chan- 
nels are changed significantly during two spring months 
(March-April) over the Kyiv site, which is confirmed by 
the Angstrom exponent variation from 1.35 to 1.50 (see in 
Fig. 2d). These spring changes are not so prominent over 
the Minsk site. In contrast the aerosol particle spectral char- 
acteristics’ variation during summer months (July- August 
over Minsk and June-August over Kyiv) is even stronger: 
the Angstrom exponent mean value difference equals about 
0.25 (over Minsk) and 0.17 (over Kyiv), which is evidence 
of a less homogeneous composition of aerosol particles in 
the summer than in the spring. The POLDER AOD (865 nm) 
values for the Kyiv and Minsk areas are shown in Fig. 2a and 
c by open diamonds, and they are lower than the AERONET 
data, having traits in general similar to the AOD (870 nm) 
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(a) (b) (a) (b) 



Figure 2. (a-c) The AOD (440 nm) and the AOD (870 nm) mean 
value seasonal changes by AERONET data and the AOD (865 nm) 
by POLDER data for the Kyiv and Minsk sites, (d) the set of 
Angstrom exponent (440-870 nm) values, (e) the features of AOD 
(440 nm), and (f) Angstrom exponent (440-870 nm) season varia- 
tions at the Kyiv site from year to year. Vertical bars show monthly 
mean standard deviations from averaged values. Months are shown 
by their first letters in the order January-December. 


values. The AOD (865 nm) is lower than the AOD (870 nm) 
because POLDER data are determined by fine-mode parti- 
cles in contrast to AERONET total (fine and coarse modes) 
particle AOD. Results in Fig. 2a-c show the season behav- 
ior of AOD and Angstrom exponent parameters on average 
having the strong variation from year to year (see Fig. 2e 
and f). In the case of the Kyiv site, AOD (440 nm) and 
Angstrom exponent seasonal variability are influenced by 
several strong events, which are not always similar in the 
AOD and the AE. For example, the strong AOD (440 nm) 
event in April 2009 corresponds to the high AE values, but 
a similar AOD (440 nm) event in April 2011 and the strong 
AOD (440 nm) event in August 2010 correspond to the mean 
AE values. 

The monthly mean total, fine-mode and coarse-mode vol- 
ume concentration values retrieved from the AERONET 
database and averaged over 2008-2013 have been used 
for estimation of aerosol particle content seasonal changes 
(Fig. 3). We could expect strong correlation between the 
total volume concentration and the AOD values. However, 


Figure 3. Seasonal behavior of aerosol volume concentration of to- 
tal, fine- and coarse-mode particles over the (a) Kyiv and (b) Minsk 
AERONET sites. 


the correspondence of these values exists only in general in 
the case of the Kyiv site. The AOD and the total volume 
concentration over the July-November period reaches maxi- 
mum in August and decreases monotonously in September- 
November. The total volume concentration over the Kyiv 
area has its minimum in winter months when fine-mode par- 
ticles prevail and increases during spring (from March to 
June) mainly due to increase of the coarse-mode particle 
concentration. During May and June the coarse-mode par- 
ticles prevail over the Kyiv site, but the AOD variation in 
both these months and in the rest of the year corresponds 
mainly to the fine-mode particle volume concentration. The 
AOD (440 nm) and AE (440-870 nm) have the sharp max- 
ima in April and August, whereas volume concentrations of 
the fine- and coarse-mode particles are practically the same 
at this time, which can be interpreted as a prevailing influ- 
ence of the fine-mode particles on the optical properties of 
the aerosol layer over the Kyiv site. 

The AOD and the total volume concentration of the parti- 
cles over the Minsk site better fit one another (Figs. 2c and 
3b), as well as in the case of fine- and coarse-mode variations. 
The relevant data have the same trend during all months ex- 
cept July and December. Total volume concentration has two 
maxima in April-May and August and has the minimum in 
winter months like in the Kyiv site. 

The difference in the AOD and the total volume con- 
centration variations over the Kyiv site can partially be ex- 
plained by conditions of the measurement selection, because 
the AOD is determined from direct sun observations that 
are performed every 15 min at SZA 28-82° during sum- 
mer (SZA >73° in winter) at the Kyiv site, but sky almu- 
cantar measurements are carried out for limited SZA values 
(Holben et al., 1998). Therefore, the volume concentration 
and physical properties of aerosol particles are determined in 
different conditions and using a smaller number of measure- 
ments. We have created the shortened series of the monthly 
mean AOD and AE measurements that consists of direct 
sun observations coincident with sky almucantar measure- 
ments to reveal the observation condition difference impact 
on seasonal variations (Fig. 2d and f). We reveal that obser- 
vational selection, i.e., amount and condition of observations, 
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influences the AOD data to a different extent for different 
months. The influence is largest in the winter and in the early 
spring. Also it is noticeable during the second half of the 
summer and in the late fall. For example, in January the av- 
eraged AOD (440 nm) value determined from the direct sun 
irradiance measurement is ~ 30 % larger than the AOD Aim 
(440 nm) value determined from direct sun irradiance mea- 
surements coincident with the sun almucantar scanning of 
the sky. Also, in March the AOD (440 nm) is ^26.5% 
larger than AOD Aim (440 nm). In summer and in two autumn 
months the AOD (440 nrn)-AODAi m (440 nm) difference is 
less than ~ 10 %, increasing to ~ 16 % in November and De- 
cember. The AOD (440 nrn)-AODAim (440 nm) relationships 
vary in the same way as AOD values increasing in the spring 
with maximum in April, and increasing in the late summer 
with maximum in August. These regularities appear at both 
the 440 and 870 nm spectral channels. The AE (440-870 nm) 
values have a reduced amount and correspond to bad obser- 
vation conditions mainly during winter. For the remainder of 
the year the AE values of both sets vary in a similar man- 
ner but with higher values of almucantar coincident AE; see 
Fig. 2d. 

Using POLDER data we investigated seasonal variation 
of the AOD (865 nm) and Angstrom exponent (670-865 nm) 
values over the Ukrainian urban-industrial areas centered on 
the cities Kyiv, Lugansk, Donetsk, Kharkiv, Dnipropetrovsk, 
Lviv, and Rivne and on the Belarus city Minsk (Fig. 4). 

Similar to the central part of Ukraine (the Kyiv site) 
and Belarus (the Minsk site) the seasonal AOD (865 nm) 
variations in the eastern (Donetsk, Lugansk, Kharkiv, 
Dnipropetrovsk) and western (Lviv, Rivne) regions of 
Ukraine exhibit two main peaks: the first in the spring 
(March-April) and the second in the late summer (July- 
August). The lowest AOD (865 nm) values are seen in June 
and in late autumn (Fig. 4). The behavior of AOD (865 nm) 
in the western part of Ukraine (Lviv) is different and ex- 
hibits earlier appearance of spring and summer AOD maxima 
(Fig. 4d). 

The comparison of the AERONET AOD (870 nm) fine- 
mode values with the POLDER AOD (865 nm) values over 
the Lugansk site shows their consistent levels and behaviors 
(Fig. 5). In the Lugansk site the two peaks in season vari- 
ation of AOD (870 nm) in total particles are determined by 
spring maxima in fine and coarse modes and by the July- 
August maximum of fine-mode (Fig. 5) by the AERONET 
data. The coarse-mode particle AOD (870 nm) observed over 
the Lugansk site has one spring season maximum in April- 
May and minimum in the late autumn (November). 

Summary statistics of the all measurements for the AOD 
and Angstrom exponent over different urban-industrial sites 
and areas of Ukraine are reported in the Supplement Ta- 
ble S 1 , which describes the averaged aerosol load. 
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Figure 4. Seasonal variations of AOD (865 nm) values and AE 
(670-865 nm) values averaged over 2008-2011 in the atmosphere 
over the Kyiv, Minsk, Donetsk, Lugansk, Dnipropetrovsk, Kharkiv, 
Lviv, and Rivne urban-industrial areas based on POLDER data. 


3.2 SSA and refractive index seasonal variability by the 
Kyiv site data 


The SSA determines the ability of the aerosol particles to 
scatter the light, and therefore it determines their influence 
on the solar energy propagation in the Earth’s atmosphere. 
SSA depends on the particle size, shape, chemical composi- 
tion, and structure. The complex refractive index is defined 
by the particle chemical composition and structure. The RI 
together with parameters of particle size distribution are the 
independent variables in the equations that describe the solar 
radiation transfer through the atmosphere, and they are deter- 
mined from the sky almucantar and direct sun AOD measure- 
ments in the sun-sky radiance inversion algorithm (Dubovik 
and King, 2000; Dubovik et al., 2006). The SSA is calcu- 
lated using the aerosol absorption and scattering coefficients 
determined from the RI and size distribution. Therefore, the 
SSA and RI parameters are physically and algorithmically 
dependent variables. 

The estimation of the SSA and RI variations by the Kyiv 
site data presents the first results for evaluation of aerosol 
spectral characteristics in the region. The level 2.0 data, 
which are appropriate for SSA and RI calculation, were not 
sufficiently collected at the Kyiv site because of weather con- 
ditions during period of the observations. For example, the 
SSA and RI level 2.0 data are absent at the Kyiv site for 
November-January, one value of the SSA and one of the RI 
have been obtained only in February, and two values have 
been obtained only in October during five years of observa- 
tions. In contrast, from tens to hundreds of the SSA and RI 
level 1.5 values are available for each year. Therefore, we 
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Figure 5. The AOD season variation over Lugansk by the AOD 
(870 nm) for total, fine and coarse mode from the AERONET site 
data and the POLDER AOD (865 nm) data. 


have to use level 1.5 for analysis of the SSA and RI seasonal 
variability. 

We averaged SSA and RI by months over the 2008-2013 
period, and it could be expected that possible errors due to a 
drift of the sun-photometer characteristics would be consid- 
ered as random values taking into account the replacement 
of sun photometers at the site each year. Therefore we as- 
sume that level 1.5 data are more or less acceptable to reveal 
seasonality of the SSA and the RI. To estimate reliability of 
that assumption, we calculated standard deviation (SD) of 
SSA values and have determined that SD of SSA (440 nm) 
is ±0.09 (11%), SD of SSA (1020nm) is ±0.12 (15%) in 
February, and SD is ±0.02 (2.5 %) and ±0.05 (6.5 %) in July 
and August, respectively. Resting SD values of the SSA are 
between these thresholds levels. 

The SSA seasonal peculiarities are seen as the spectrum 
difference between the warm period (April-October) and 
cold period (November-March). In the warm period, the 
SSA (440 nm) value is within the range of 0.89-0.93, SSA 
(1020 nm) within 0.84-0.89 (Fig. 6a and b), and the SSA 
decreasing versus wavelength is less sharp than in the cold 
period. In the cold period the SSA (440 nm) value is within 
0.86-0.82, and the SSA (1020 nm) values within 0.79-0.75, 
with SD values of ±0.07 (8 %) to ± 0.10 (13 %). 

To evaluate the effect of observational selectivity we cal- 
culate monthly averaged SSA values for the spring-summer 
period using level 2.0 data (Fig. 6b). These SSA values are 
less homogeneous for different wavelengths than level 1.5 
data. However, the SSA level 1.5-level 2.0 value difference 
is within the SD range of the SSA level 1.5 except for the 
March data, which are explained by the difference in the 
amount of observations. 

Similarly to the SSA parameter behavior, the RI imagi- 
nary part (Rlim) spectrum exhibits noticeable seasonal dif- 
ference between warm and cold parts of a year (Fig. 6c). 
The highest monthly averaged RIi m values are0.018±0.017 
at 440 nm and 0.022 ±0.024 at 1020nm in the warm pe- 
riod. In June, RIi m values are the lowest: 0.009 ± 0.005 at 
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Figure 6. (a and b) Single-scattering albedo, and (c and d) refractive 
index spectral dependence by the Kyiv AERONET site data. 


440 nm and 0.010 ±0.008 at 1020 nm. During November- 
March, aerosol particles are defined by relatively large RIi m 
values: 0.027-0.037 at 440 nm and 0.034-0.051 at 1020 nm. 
Since RIi m defines the absorption ability of the aerosol par- 
ticles, the RIi m spectrum seasonal behavior is antisymmetric 
to the SSA spectrum dependence (see Fig. 6a and c). 

The RI real part (RI re ) in the average shows seasonal vari- 
ations similar to the RIi m values but less distinct (Fig. 6d). 
The RI re monthly averaged spectral characteristics’ values 
during December-March are the lowest in the whole RI re 
range. The shapes of the averaged RI re spectral character- 
istics during this season are similar to the RIi m spectrum and 
are mirror-like to the SSA spectrum behavior. This tendency 
is typical for urban-industrial aerosols according to Dubovik 
et al. (2002). 


4 Discussions 

4.1 Seasonal aerosol variations over the Ukraine 
urban-industrial areas 

Several investigations of aerosol events in the eastern Euro- 
pean region reported as the case studies (e.g., Arola et al., 
2007; Zawadzka et al., 2014). In our paper we are trying to 
describe the seasonal aerosol behavior integrated over sev- 
eral years period for Ukrainian regions that is expected to be 
exposed by urban-industrial aerosol load. 

The first spring AOD peak seen in the almost all areas by 
POLDER data partly can be explained by natural changes 
in land surface between snow cover in winter and open soil 
in other seasons, which is accompanied by an increase of 
soil dust (Sterk and Goossens, 2007; Hinz and Funk, 2007), 
traffic pollution, and biomass burning aerosol (Stohl et al., 
2007; Bamaba et al., 2011) in the atmosphere. The second 
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(July-August) peak in AOD values appeared in the averaged 
data due to frequent wildfires that occurred in the years con- 
sidered (2008, 2010, 2011; see Fig. 2e). The summer AOD 
peak for the western part of Ukraine (see Fig. 4d) is shifted to 
the midsummer (July), relatively to those in eastern and cen- 
tral Ukraine, which agrees well with the results reported by 
Jaroslawski and Pietruczuk (2010) for the Belsk site (Poland) 
with fine-mode particle domination (Dubovik et al., 2011). 
This summer maximum is related to the seasonal biomass 
burning started earlier in eastern and southern Europe. 

According to the AOD (865 nm) data from the POLDER 
instrument, the higher levels of aerosol contamination are ob- 
served over the Kyiv and Lviv cities in the winter/spring and 
over Kharkiv in the summer. Surprisingly twice-lower AOD 
(865 nm) values were observed over industrial areas Donetsk 
and Lugansk. However it could be partly explained by se- 
lectivity of POLDER data that are mainly represented by 
fine-mode aerosol. Contamination by coarse-mode fraction 
in these areas can be significant. This result was confirmed by 
the Lugansk AERONET site measurements (Fig. 5), where 
the total aerosol AOD (870 nm) value is almost the same as in 
the Kyiv site. The AERONET measurements at the Donetsk 
site, made episodically in the summer, confirm the significant 
coarse aerosol mode presence in the data. 

The Angstrom exponent (670-865 nm) averaged values 
have a seasonal trait with two maxima, the first in the spring 
(April) and the second in the summer (July-August), and a 
minimum in late autumn for Kyiv, Donetsk, Dnipropetrovsk, 
Kharkiv, and Lviv areas (see Figs. 2d and 4). In Minsk, Lu- 
gansk, and Rivne areas the AE values vary in different ways, 
increasing from winter months to maximum in midsummer 
and decreasing in autumn. Therefore we can expect that on 
average the fine-mode particles are dominant in these areas 
in summer months. The comparison of AOD (870/865 nm) 
for the Kyiv and Minsk AERONET sites’ data and POLDER 
data for these two areas shows general agreement in vari- 
ations and values. There is a difference between AE (440- 
870 nm, AERONET) and AE (670-865 nm, POLDER) be- 
havior over the Minsk area (Figs. 2d and 4a) with a differ- 
ent period of the summer AE maximum appearance. The AE 
(440-870 nm) in Fig. 2d was retrieved from the AERONET 
total aerosol particles amount, while the AE (670-865 nm) in 
Fig. 4a was calculated from the fine-mode particle fraction 
observed by the POLDER instrument. The different spectral 
pair line AE ratio and the AE dependence on the bi-modal 
particle size distribution have been investigated by Schus- 
ter et al. (2006). They elucidated that spectral variation of 
aerosol extinction is important in calculations of the AE with 
two wavelengths only. The results of the AE calculation from 
the longer wavelength pair (670-870 nm) are sensitive to the 
fine-mode aerosol fraction and do not depend on the fine- 
mode particle effective radius. The AE calculated from the 
shorter wavelength pair (380-440 nm) is more sensitive to 
the fine-mode effective radius than to the fine-mode fraction 
(Schuster et al., 2006). Therefore, the AE (440-870 nm)-AE 


(670-865 nm) relationship is explained by the influence of 
different particle properties. 

The amplitude of the seasonal variations in the Sevastopol 
AERONET site is much smaller in comparison with the Kyiv 
and Minsk AERONET sites and other investigated Ukrainian 
areas (Bovchaliuk et al., 2013). This feature is explained by 
coastal position of the Sevastopol site at the Black Sea and 
less strong urban and industrial impact than in other Ukraine 
regions. 

4.2 Seasonal aerosol behavior over the Kyiv site 

The data analysis shows that the aerosol AOD values and 
volume concentrations in the atmosphere over the Kyiv site 
are on average larger than over the Minsk site, whereas the 
AOD seasonal variations are similar (Figs. 2-4). The AOD 
over both sites exhibits a prominent peak in August and sim- 
ilar variations over the autumn and the winter, with signif- 
icant deviations from the annual mean in individual years. 
This seasonal variation is evidence of similarity of the basic 
aerosol sources in these two regions and is closely related to 
properties of aerosol particles. 

The increased AOD values in the spring and the late sum- 
mer over the Kyiv site correspond to a relative increase of the 
aerosol fine-mode particle concentration during each year. 
In contrast the coarse-mode relative content increases when 
the AOD decreases in midsummer and winter (see Fig. 3). 
The AOD (440 nm) maxima both in April and in August ap- 
pear when the volume concentration values of the fine- and 
coarse-mode particles are approximately equal (Figs. 2 and 
3), but in general the AOD (440 nm) and the AOD (870 nm) 
seasonal changes better fit to the fine-mode particle volume 
concentration seasonal behavior, which explained the SSA 
and RI values’ seasonal changes as well. The AE (440- 
870 nm) variation indicates domination of the fine-mode par- 
ticles in the AOD (440 nm) maxima, especially in those that 
are prominent in Fig. 2d, where the AE was retrieved us- 
ing the AOD almucantar measurements. The inflection point 
between fine and coarse mode on the particle volume size 
distribution (Fig. 7) is varied during the year from ^0.5 
to ~ 0.9 pm with minimum in May over the 2008-2013 pe- 
riod at the Kyiv site (not shown). In general, the fine-mode 
aerosols define the AOD behavior over the Kyiv site during 
the year except during May-June (Fig. 2a). 

The AE values over the Kyiv site are highest in the mid- 
summer, which can be explained by the domination of the 
accumulation mode particles in the AOD. The latter deter- 
mines the AE behavior (see Fig. 2a) due to summer biomass 
burning and urban-industrial aerosols increasing (Eck et al., 
1999). It should be emphasized that the amount of coarse- 
mode particles increases in the period from January to June 
(Fig. 3a). The noticeable feature is that AOD (440 nm) values 
in the beginning of the summer (May-June) have practically 
the same values for each year (Fig. 2e), but the AE values are 
different (Fig. 2f). 
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Figure 7. Monthly mean aerosol particle volume size distribution in the 2008-2013 period for the months that describe (a) maximum AOT 
(440 nm) in spring (April, filled symbols) and minimum in winter (January, open symbols); (b) maximum AOT (440 nm) in late summer 
(August, filled symbols) and minimum in autumn (November, open symbols). The Kyiv AERONET site level 2.0 data. 


The Angstrom exponent seasonal changes are similar at 
the Kyiv and Minsk sites, with a peak in midsummer (Kyiv) 
and the late summer (Minsk), and a minimum in November. 
These results correspond to fine- and coarse-mode particle 
behavior in (Fig. 3), where coarse-mode prevails in May- 
June over both sites and the fine- and coarse-mode AODs are 
similar in July-September. 

The comparison of the SSA and RI data from the Kyiv 
AERONET site with results for other European regions al- 
lows concluding that aerosol particles over the Kyiv site are 
of the urban-industrial type according to classification of 
main aerosol types (e.g., Dubovik et al., 2002; Bergstrom 
et al., 2007; Giles et al., 2012). The seasonal variations of 
the SSA and the RI represent particle composition changes; 
in particular the AOD increase in the spring due to coarse- 
mode particle concentration rise (Figs. 3 and 5) is connected 
to soil dust raising that corresponds to the raised RI re values 
in April (Fig. 6d). 

The analysis of aerosols above Moscow during wildfires 
in 2010 done by Chubarova et al. (2012) shows that smoke 
aerosols from forest wildfires also influence the SSA and RI 
values and their spectral characteristics. That sort of events 
happened as usual in the warm season. In our analysis we do 
not exclude these events in data averaging. It should be noted 
that the SSA and RI are loaded with considerable systematic 
biases because sky-radiance measurements are impossible on 
cloudy days. The additional biases arise from break inversion 
stability criteria that require high aerosol loading conditions 
when the RI re , RIi m , SSA retrieval can be more accurate. The 
aerosol content above the Kyiv site is relatively small: the 
AOD (440 nm) exceeds 0.4 only during the special events, 
like forest fires in 2010. 

The monthly mean volume size distributions for the 
months that represent the season peculiarities - maximum 
AOT (440 nm) in the spring (April) and minimum in the 
winter (January), maximum AOT (440 nm) in the late sum- 
mer (August) and minimum in the autumn (November) over 
2008-2013 - are presented in Fig. 7. The significant differ- 
ence in the coarse-mode of the particle size distribution is 
seen between January and April, and November and August. 


In the winter and autumn months the fine-mode particles with 
modal radii of about 0.2 pm prevail and coarse-mode parti- 
cles are almost not observed. In April over the Kyiv site, the 
dV (r)/dlnr value of the fine-mode particles is 3 times larger 
than in January and the particle modal radii is ~ 0. 1 3 pm. The 
coarse-mode particles (2-5 pm) with similar content are ob- 
served in April and August as well. 

To find the potential types of aerosol over the Kyiv 
AERONET site we created scatter plots (see Fig. 8) of 
the SSA versus the extinction Angstrom exponent (EAE) 
for April and August months separately, the absorption 
Angstrom exponent versus the EAE (440-870 nm), the 
Angstrom exponent (440-870 nm) versus the AOD (440 nm), 
and the Angstrom exponent difference (AE (440-675 nm)- 
AE (675-870 nm)) versus AE (440-870 nm) according to 
classifications developed in Dubovik et al. (2002), Omar et 
al. (2005), Russel et al. (2010), Giles et al. (2012), Boselli et 
al. (2012), and Yoon et al. (2012). 

Analysis of scatter plots in Fig. 8a and b shows that aerosol 
arrived/measured over the Kyiv AERONET site belongs to 
the mixture of urban-industrial type mostly and partly to the 
biomass burning type. Dust particles are not identified ac- 
cording the clustering by aerosol type given in Giles et al. 
(2012). The seasonal discrimination of Angstrom exponent 
difference (AE (440-675 nm)-AE (675-870 nm)) versus AE 
(440-870 nm) scatter plot may allow dividing particles by the 
type prevailing in each season (Fig. 8c and the Supplement 
Fig. SI). 

According to a classification proposed by Yoon et 
al. (2012) the fine-mode urban-industrial aerosols in win- 
ter months (DJF) are observed with small presence of dust 
coarse-mode particles. In spring months (MAM) the urban- 
industrial type and dust coarse-mode aerosol are observed in- 
cluding a smaller-amount biomass burning type. In summer 
months (JJA) the fine-mode urban-industrial aerosol is dom- 
inant and the biomass burning type aerosol is observed in 
contrast to SON. In autumn months (SON) urban-industrial 
and mixed type of aerosol is dominant; however, the other 
type of aerosol particles appeared in the SON period only 
in segment with the AE-difference / AE ratio equal to about 
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Figure 8. (a) Scatter plot SSA versus extinction Angstrom expo- 
nent (EAE) for April and August months separately, (b) absorp- 
tion Angstrom exponent versus EAE (440-870 nm). (c) Seasonal 
scatter plot Angstrom exponent difference (AE (440-675 nm)- 
AE (675-870 nm)) versus Angstrom exponent (440-870 nm), and 
(d) Angstrom exponent (440-870 nm) versus AOD (440 nm). Data 
from the Kyiv AERONET site. 


1.5 / 1.5, which is not included in the Yoon et al. (2012) clas- 
sification (see Fig. 8c and Fig. SI). 

The analysis of the seasonal differentiation of the AE 
(440-870 nm) versus AOD (440 nm) scatter plot allows as- 
suming the seasonal features of aerosol particles over Kyiv 
(Fig. 8d and the Supplement Fig. S2). In DJF the AOD 
(440 nm) values are on average less than ~0.4. The fine- 
mode particles with AE> 1.8 are not observed in the win- 
ter season. In the spring only the coarse-mode particles with 
AE<0.5 and AOD (440 nm) ~ 0.1 -1.1 are measured, and 
fine-mode particles with AE ~ 1.8 and AOD (440 nm) up to 
0.9 are measured. A small amount of the fine-mode frac- 
tion dust with AE>2 is observed in summer months. The 
event of wildfire in Russia (Chubarova et al., 2012) in mid- 
August 2010 seen well as a chain of data with AE~ 1.4 
and AOD (440 nm) up to ~ 1.5 as the result of the biomass 
burning smoke particles arrived at the Kyiv site on 15- 
16 August 2010. The autumn AE/AOT (440 nm) distribu- 
tion is similar to the spring scatter plot except for the smaller 
amount of fine-mode (AE> 1.8) and coarse-mode (AE<0.5) 
particles. The scatter plots data correspond to results of the 
AOD and Angstrom exponent seasonal traits in Fig. 2: min- 
imal AOD values in the winter and in the autumn and maxi- 
mum ones in the summer. 

We applied the back-trajectory and cluster analysis meth- 
ods to calculate the seasonal back trajectories and prevail- 
ing directions of the arrived air mass for the Kyiv and Minsk 
sites to investigate the possible seasonal variations of aerosol. 


The 3-D back trajectories over the 2012-2013 period for the 
Kyiv site and over the 2008-2013 period for the Minsk site 
have been downloaded from the AERONET website. The 
back trajectories over the 2008-201 1 period for the Kyiv site 
we calculated using the software developed in Burak et al. 
(2010). The wind field and temperature data from the Be- 
larus Hydrometeorology Center (http://hmc.by/) were used 
for calculation. 

The calculated trajectory presents the pathway of air mass 
transported at the height of about 1.5 km (850 hPa) dur- 
ing 48 h before measurement. The cluster analysis produced 
the five cluster trajectories with prevailing directions aver- 
aged by seasons over the period 2008-2013. The trajecto- 
ries in day of measurement (red lines) at the Kyiv and Minsk 
AERONET sites and the seasonal cluster trajectories from 
five prevailing directions (bold blue lines) are shown in the 
Supplement Fig. S3. For each season we determined the most 
prevailing trajectory that is responsible for arrived aerosol. 
This trajectory is described by fraction N of the trajectories 
number in the main cluster to the number of trajectories in 
all five clusters. The most prevailing trajectories for the Kyiv 
and Minsk sites in each season are shown in Table 2. 

We consider a possibility to determine the presence of the 
Saharan dust by the Kyiv and Minsk site data. According 
to Israelevich et al. (2012) the events with the high AOD 
values in the European region are associated with biomass 
burning, anthropogenic pollution, soil erosion, and trans- 
ported dust from Saharan sources. All of these events are ex- 
pected to exhibit seasonal variations. For the eastern Mediter- 
ranean three periods of increased atmospheric dust were ob- 
served: in spring (March-May), in summer (July-August), 
and in autumn (September-November). Results of ten years 
(2001-2010) of MODIS data analysis of the seasonal vari- 
ations of Saharan dust transport over the Mediterranean to 
Europe (Israelevich et al., 2012) show AOD seasonal in- 
crease in April and July-September at the Kyiv longitude 
similar to our results in Fig. 2. Papayannis et al. (2008) re- 
ported on Saharan dust event observations on about 20 days 
by Minsk lidar measurements over the 2000-2002 period. 
Saharan dust intrusions in the eastern Mediterranean and 
Europe are most active in the spring season (Barnaba and 
Gobbi, 2004; Papayannis et al., 2008). Therefore we can ex- 
pect the appearance of Saharan dust to the Kyiv and Minsk 
sites in 2008-2013. According to seasonal back trajectories 
the most favorable conditions for dust aerosol transport from 
the eastern Mediterranean region exist in the spring and in 
the autumn. The spectral aerosol parameters SSA and RI in 
spring months (Fig. 6) are close to those reported by Kim 
et al. (2011) for dust aerosol. Moreover, the SSA spectrum 
in May 2008 (not shown) exhibits similar spectral variation 
and values to the desert dust aerosols classified by Dubovik 
et al. (2002). Aerosol particles with parameters AE<0.5 
and AOD (440 nm) >0.4 are observed at the Kyiv site in 
the spring and summer months (Fig. 8d and the Supple- 
ment Fig. S2b and c) that correspond to dust type aerosols 
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Table 2. The most prevailing seasonal trajectory for the Kyiv and Minsk sites calculated using the back trajectories and cluster analysis. 


Season/ AERONET site 


Kyiv 

Minsk 

Direction 

Fraction N 

Direction 

Fraction N 

Winter 

North 

25/87 

Northeast 

31/107 

Spring 

West 

223/753 

Northwest 

158/570 

Summer 

West 

170/709 

West 

125/482 

Autumn 

South 

167/424 

East 

57/159 


described by Pavese et al. (2012). Therefore the Saharan min- 
eral dust that often extended across the Mediterranean into 
Europe (Papayannis et al., 2008; Israelevich et al., 2012) is 
transported to the eastern European region mostly in the late 
summer and autumn (August-September) and infrequently 
in the spring (April). This inference is confirmed by increas- 
ing of the RI real part in these months that corresponds to the 
mineral particle composition increase (Fig. 6). 


5 Conclusions 

The seasonal variations of aerosol amount and optical prop- 
erties over several Ukrainian urban-industrial areas are ana- 
lyzed using the ground-based data from AERONET sun pho- 
tometers and POLDER satellite instrument data. We investi- 
gate the aerosol seasonal behavior over the areas of the cities 
Kyiv, Lugansk, Donetsk, Kharkiv, Dnipropetrovsk, Lviv, and 
Rivne and the area of the Belarus city Minsk for compar- 
ison to Kyiv. These areas are influenced by numerous lo- 
cal aerosol pollution sources (intensive city transport, heavy 
industry) and by the aerosol transport from remote sources 
(open steppe fields, mining, and wildfires). 

We extracted the monthly averaged AOD at 440 and 
870 nm and the Angstrom exponent values for the 440- 
675 nm and 440-870 nm channels from the AERONET 
database for the Kyiv, Minsk, and Lugansk sites. Also the 
size distribution and volume concentration, the SSA and 
complex RI, the absorption AOD and the extinction AOD 
retrieved by AERONET inversion algorithm were used for 
aerosol behavior analysis at the Kyiv site. The monthly mean 
AOD values at 865 nm and the AE values at 670/865 nm were 
derived from POLDER measurements. The comparison of 
the AERONET and POLDER data shows that AOD (865 nm) 
determined from POLDER measurements are lower than 
AOD (870 nm) from AERONET measurements because 
POLDER data are determined by fine-mode particles in con- 
trast to AERONET total (fine and coarse modes) particle 
AOD. 

We reveal the similar behavior of average aerosol param- 
eter values in all areas of the eastern and central part of 
Ukraine where measurements were performed using both 
AERONET and POLDER techniques: the aerosol amount 
and the optical thickness values exhibit peaks in the spring 


(April-May) and the late summer (August), whereas mini- 
mum values are seen in the late autumn. The summer AOD 
maximum is larger in the Kyiv and Minsk sites, and min- 
imum aerosol content is reached in the late autumn. The 
behavior of the aerosol parameters in the western part of 
Ukraine is different from eastern and central regions and ex- 
hibits an earlier appearance of spring and summer AOD max- 
ima. This aerosol behavior in general agrees with the results 
achieved by other authors, e.g., Sayer et al. (2012) on sea- 
sonal variations in the Eurasian region and by Pietruczuk and 
Chaikovsky (2012), who suggested also that biomass burn- 
ing products and Saharan dust are responsible for large AOD 
values during the spring time in this area. 

The aerosol particle volume concentrations over the Kyiv 
and Minsk sites have two maxima during the warm part of 
a year. The AOD and the total, as well as fine- and coarse- 
mode volume concentration of the particles over the Minsk 
site have the same as in the Kyiv site trend during all months 
except July and December, and two maxima: in April-May 
and in August. The correspondence of these values exists in 
general also in the case of the Kyiv site, but the total particle 
volume concentration over Kyiv reaches maximum in June 
simultaneously with the volume concentration of the coarse- 
mode particle maximum, whereas the AOD maximum comes 
in April simultaneously with the fine-mode particle volume 
concentration maximum. The total volume concentration has 
its minimum in winter months, when fine-mode particles pre- 
vail. 

On average, the aerosol amount over the Kyiv site is larger 
than that over the Minsk site in the summer and in the early 
autumn. The load of fine-mode particle content over Minsk 
is larger than over the Kyiv site during almost all seasons. 

To reveal the aerosol particle microphysics seasonal be- 
havior and to estimate the aerosol type, the SSA and RI 
values and spectra at the Kyiv AERONET site were an- 
alyzed. The SSA monthly averaged spectra are divided 
into two groups: during warm period (April-October) they 
have higher values and relatively smooth and monotonous 
shapes, whereas during cold period the SSA values are 
lower and shapes are more tilted, bended and more vari- 
ous. The monthly averaged RI spectra also are different in 
warm and cold parts of a year. Specifically the RI imagi- 
nary part values and spectra have sharp difference for the sea- 
sons and are mirror-like to the SSA spectrum behavior. The 
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monthly averaged spectral characteristics of the RI re during 
December-March are in the lowest part of the range. The 
SSA and refractive index values and spectra correspond to a 
mix of urban-industrial, biomass burning, and dust aerosols. 

The amount of coarse-mode particles in the atmosphere 
above the investigated Ukrainian areas during the warm pe- 
riod of a year exceeds, on average, twice their amount during 
the cold period. We consider that as an effect of large open 
land areas in the south of Ukraine and in the southeastern Eu- 
ropean part of Russia where dust storms are frequent during 
the warm period, whereas during the cold period these areas 
are mostly covered by snow preventing dust uplift. 

We create scatter plots of the different aerosol parameters 
using the Kyiv site data to evaluate the potential types of 
aerosol. According to a classification proposed by Yoon et 
al. (2012), in winter months the fine-mode urban-industrial 
aerosols are observed with small presence of the dust coarse- 
mode particles. In spring months the urban-industrial type 
and dust coarse-mode aerosol are observed including the 
smaller-amount biomass burning type. In summer months 
the fine-mode urban-industrial aerosol is dominant and the 
biomass burning type aerosol is observed. In autumn months 
urban-industrial and mixed type of aerosol is dominant; 
however, the other type of aerosol particles appeared in the 
SON period only in segment with the AE-difference / AE ra- 
tio, which is not described in the Yoon et al. (2012) classifica- 
tion. The scatter plot data in general correspond with results 
of the AOD and Angstrom exponent seasonal traits, with the 
low AOD values in the winter and in the autumn and large 
AOD in the spring and the summer. This seasonal behavior is 
explained by the influence of the local sources as well as by 
aerosol transport. According to back-trajectory analysis the 
industrial type aerosol particles are transported mainly from 
the western Europe regions in the spring and in the summer. 

The presence of the Saharan dust was investigated using 
the Kyiv site data. According to seasonal back trajectories the 
most favorable conditions for dust aerosol transport from the 
eastern Mediterranean region, which experienced frequent 
Saharan dust intrusions, exist in the spring and the autumn. 
The spectral aerosol parameters (SSA, RI) measured at the 
Kyiv site in spring and in summer months comply with the 
dust type particles (Kim et al., 2011; Pavese et al., 2012). 
Therefore, we can assume that the Saharan mineral dust ex- 
tending across the Mediterranean into Europe (Israelevich 
et al., 2012) is possibly transported to the Ukrainian areas 
mostly in the late summer and autumn (August-September) 
and infrequently in the spring. 


The Supplement related to this article is available online 
at doi:10.5194/amt-7-1459-2014-supplement. 
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